Efficient excitation of a two level atom by a single photon in a propagating mode 
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State mapping between atoms and photons, and photon-photon interactions play an important role in scalable 
quantum information processing. We consider the interaction of a two-level atom with a quantized propagating 
pulse in free space and study the probability P e (t) of finding the atom in the excited state at any time t. This 
probability is expected to depend on (i) the quantum state of the pulse field and (ii) the overlap between the 
pulse and the dipole pattern of the atomic spontaneous emission. We show that the second effect is captured by 
a single parameter A 6 [0, 8w/3], obtained by weighting the dipole pattern with the numerical aperture. Then 
P e (t) can be obtained by solving time-dependent Heisenberg-Langevin equations. We provide detailed solutions 
for both single photon Fock state and coherent states and for various temporal shapes of the pulses. 
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I. INTRODUCTION 

Light-matter interface in free space at quantum level lies at 
the heart of quantum networks and quantum communication 
as well as being the fundamental question in quantum 
optics, and may be less technologically demanding than 
typical cavity quantum electro-dynamical systems. Recently, 
achievements have been made for atoms fl-@]> single 
molecules |0, 01 and quantum dots Hfl. Especially, high 
numerical aperture optics [9] have been recognized as a key 
element for free space atom-light coupling andprecision 
spectroscopy, involving fixed aspheric lenses @, Bl fT^ - H^ll . 
deep parabolic mirrors lfl3T[l5ll . spherical mirrors lUotl and 
phase Fresnel lenses lfl7ll . 

Regarding the atom-light interaction, there are two different 
phenomena: the scattering of the light by the ato m pi EI [ill an d 
the absorption of the light by the atom H, [H, Efficient 
information transfer between atoms and photons requires 
controlled photon absorption with high probability. However, 
the time reversal argument implies the main properties of the 
excitation pulse. These are (i) the spatial profile of the pulse 
should match the atomic dipole emission pattern, and (ii) the 
temporal shape of the photon has to be a time inverted version 
of the spontaneously emitted photon. This means that the 
atom must be illuminated from all directions by a photon 
with a rising exponential temporal envelope fl5l[l8ll . Strong 
focusing can give rise to increased overlap of the light beam 
with the atomic dipole and this improvement of the atom-light 
coupling in free space has been predicted by theory I il l 1(111 and 
demonstrated in experiment 141 la [Till . 

In this paper, we focus on the effect of the temporal-spectral 
features of the pulse on the probability of finding the atom in 
the excited state starting from the ground state ("excitation 
probability"). We present the general formalism and apply it 
into two specific situations: single photon wave packets and 
coherent state wave packets. 

The paper is organized as follows. In Sec. [TTJ we review 
a general quantized model of the interaction between an atom 
and a continuum propagating pulse in free space and introduce 
the basic parameters describing the coupling between the 
atomic dipole and the pulse. In Sec. [Till we introduce a 
special case. In Sec. [IV] the dependence of the atomic 
excitation probability on the temporal and spectral features of 



single photon wave packets and coherent state wave packets 
are investigated, respectively. The excitation probability for 
realistic focusing setups is discussed in Sec[V] And our results 
are briefly summarized in SeclVIl 



II. GENERAL MODEL AND APPROACH 

We start by considering the interaction of a two-level atom 
sitting at position r a with the quantized radiation field with 
continuum modes in free space. In Coulomb gauge, the 
positive-frequency parts of the electric field operators can be 
expanded as 11911 . 
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where = c|k|, c is the vacuum speed of light, eo is the 
permittivity of the vacuum, e^iA = 1,2) are unit polarization 
vectors, ■ = 6 AA >, ^k,i • k = 0, and the field operators 
follow the usual commutation relation 



[a M ,fl^ = 5(k-k')d A 



(2) 



Energy conservation implies the normalization of the 
spatial mode functions u^ A {r) 
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In the interaction picture and rotating-wave approximation, 
the dynamics of the system is described by the Hamiltonian 

H/ = -01 J] f ^ 3 k [g k Jr a ) & + a KA e- i ^-^' - h.c] , (4) 

A 

where (±> a — E e - E g is the atomic transition frequency and 
<3"+ = \e){g\,a-- = \g){e\,& z = \e){e\ - \g){g\ = &+&- - <x_<x + 
are atomic operators. The coupling strength is given by 



■ «M(r a ) erf • e k i, 
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where d is the value of the dipole momentum and is the unit 
dipole vector. 



2 



The evolution of the operators is governed by a set of 
coupled Heisenberg equations 

*M = gt A {r„)&-e iin) ', (6) 
<5- z = -T'(fr z + 1) 
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in which the decay term proportional to P and the noise 
operators £ are introduced to account for the interaction of 
the atom with the environment. 

By integrating Eq.©, the field operator is decomposed into 
a free field part and a part radiated by the atom j20l p. 393] 
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The substitution of Eq.© back into Eq.© and ([8]l gives a set 
of modified optical Bloch equations B21I1 . 

d- z = -T(fr z + 1) 

_2 Z | [^M( r fl)0"+«k„ 



(toY 



■i(a) k -a) a )t 



h.c] + t z , 



~2 a - 
+ cr 7 



<^kg M (r fl )a M (fo)e 



(10) 



(11) 



where the standard spontaneous decay rate in free space is 
made up of two parts ll22ll : T = P + T p , the decay into the 
environment T', which is the non-pulse mode in our case, and 
the decay to the pulse mode T p . 

Following the same reasoning for the environment field 
operators by^ as for the field operators a^j, one readily finds 
the explicit form of the noise operators 

i- = ^Zj ^ 3 k'«v^( r «)^k'-i'( < b)«" ,(flw " wjf . ( 13 > 

where gk\A'( r a) is the corresponding coupling strength to the 
atom. 

Furthermore, with the use of the Weisskopf-Wigner theory 
I23I p. 207], the explicit formula for the T p can be found and 
is given by, 
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Substituting for from Eq.© and going to the spherical 
coordinates, one gets 
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where the integration runs over the solid angle d£2 covered by 
the pulse mode and the indices k a respect the condition |k a | = 
k a coming from the delta distribution in Eq.(fT4b. Note, that 
in the special case, where uy^^ = e ,k °' r ° and the integration is 
performed over the whole solid angle, the parameter lambda 
reaches its maximum value A = 87r/3 and one obtains the well 
known formula for the spontaneous decay of a single atom in 
a free space (see e.g. II2(X p. 530]) 
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This is as well the maximum possible value (in principle) for 
the pulse mode decay T p = T^^. Thus the single parameter 
A e [0, 87r/3] describes which part of space, weighted by the 
atomic dipole moment, is covered by the pulse. 

In the following, we study the excitation probability P e of 
the atom excited by the photon wave-packet, which is given 
by the expectation value of atomic operator fr z 



Pe(t) = Uwomwo) + 1), 



(17) 



where the initial state of the system |^o) = |^)|O p )|0 e ) is a 
product state of the atomic ground state, the pulse state, and 
the environment. We assume the environment to be initially 
in the vacuum state. In order to find the value of (o" z (0) at 
arbitrary time t, one has to solve the set of coupled differential 
equations for all involved time-dependent operators ( in our 
case fr z and <r ± , where the equation for & + is just a hermitian 
conjugate of Eq.^TTTi ). The complete set of equations is 
obtained by first averaging the Eq. (TlOb over the initial state 
l^o). Next, the initial state dependent average of operators &± 
has to be found and can be obtained from Eqs. (fTTt and H.c. 
of (lilt . The complete set of equations can be schematically 
written as 



s(f) = Ms(f) +b. 



(18) 



The form of the vectors s, b and the matrix M is 
initial state dependent and will be specified in Sec. [Ill] 
Note that the Langevin-type noise operators Eq. ( 1121131 ) are 
determined directly in terms of the initial field operators of the 
environment II23[ p. 273], so their average values will vanish 
as < \ >=0, when considering the initial vacuum state of the 
environment. 



DDL SPECIAL CASE: DIPOLE PATTERN 

With the help of the presented general model, one can study 
the dependence of excitation probabilities on both the spatial 
and temporal properties of the pulse. In the following, we 
assume that the light field in question matches the atomic 
dipole field pattern. Therefore, all the functions gk,.i will be 
defined with the dipole pattern. Then we only focus on the 
temporal and spectral effects of pulse. We first examine the 
interaction of the atom with a single photon Fock state pulse 
and next with the coherent state pulse. 
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A. Single photon wave-packet 

Let's first consider the pulse mode to be a single photon 
wave-packet, which can be written as J23l p. 208] 

|ip> = 2 f AriU(O/fa*)akl0> 

where the spectral distribution function f(u>k) is the only 
degree of freedom left. The normalization of the single photon 
wave-packet Eq.([T9ll implies, 

2 f rf 3 k|^ M (r a )| 2 |/(^)l 2 = 1. (20) 

The explicit form of Eq.(fT8l) can now be found with the initial 
state PP > = ls>|lp>|0,>: 



definition of continuous mode coherent state presented in 
12411 . we define the coherent state as 



\a p ) = exp [fffl[. - a*a c ] |0> , 



(25) 



where the wave-packet operator a c is defined in Eq.(|T9T>. The 
mean photon number N in the wave packet is given by, 

N = (a p \ a\ac \a p ) = \a\ 2 . (26) 

This yields the usual relation for coherent state 

a c \a p ) = a \a p ) . (27) 

Again, we get a set of similar differential equations with the 
same initial condition, but different variables 
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with initial condition 



s r (f ) = (-l 0). 



Once again with the help of Weisskopf-Wigner 
approximation, where we assume that the coupling g^ 
is constant for frequencies of interest centered around the 
atomic transition frequency (±> a , one finds that the effective 
coupling strength g(t) is the product of temporal envelope of 
the pulse and the decay rate into the pulse mode, 



g (f) = fi~ P m 



where 



m = 



2k 



J daj k f(a> k )e-* n)t , 



(22) 



(23) 



is up to a constant factor the Fourier transform of the spectral 
distribution function. The excitation probability is then given 
by the first component of the vector s 

Peit) = ^l(0 + 1). (24) 



B. Coherent state wave-packet 

Let us now have a look at the wave packet initially prepared 
in a (continuous field) coherent state. In analogy with the 



(28) 



IV. ANALYSIS OF THE TEMPORAL ENVELOPE 

In this section, we assume that the incoming pulse not only 
matches the atomic dipole pattern but also occupies the whole 
solid angle, which implies that T p = T. We then discuss the 
pulse temporal shape effect on the excitation probability. 



A. Pulse bandwidth effects 

For a fixed pulse envelope, the excitation probability 
depends on ratio between the pulse bandwidth Q. and the 
decay rate T of the atomic dipole. We take a single photon 
Fock state pulse with a Gaussian temporal shape as an 
example, and study the effects of different bandwidths on the 
excitation probability. The results are plotted in Fig[T] 

As we can see from Fig[T] for single photon excitation 
with shorter pulses (Q » T), the bandwidth is too broad 
for resonant absorption, which reduces the effective coupling 
strength. For longer pulses (Q <sc F), the photon density is too 
low for efficient interactions lf2lll . 

In Fig |2] we show the dependence of the maximum 
achievable resonant excitation probability on the pulse 
bandwidth for the single photon Fock state pulse and single 
photon coherent state pulse, where the mean photon number 
equals 1. We find out that the optimum pulse bandwidth 
maximizing the absorption is Qo = 1.5F for single photon 
Fock state pulse and Q.' () = 2.4F for single photon coherent 
state pulse. 

For coherent state pulses, we studied the maximal 
excitation probability as a function of the mean number of 
photons N for various choices of the bandwidth, shown in 
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FIG. 1: (Color online) Excitation probability P e (t) as a function of 
time with the initial Gaussian pulse in single photon Fock state for 
different bandwidths. fl = 1.5r, which turns out to be the optimized 
bandwidth (c.f. FigO- 
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FIG. 2: (Color online) Dependence of maximum excitation 
probability P^ lm on the pulse bandwidth with Gaussian shape for 
single photon Fock state pulse and single photon coherent state pulse. 



Figj3] As expected, the maximal excitation probability varies 
with N. The saturation for large N for all bandwidths is due to 
the fact that the effective coupling strength g( f) decreases with 
the pulse length. Alternatively, this can be understood as the 
photons arrive more distributed in time. Note that, for large 
N, it is better to choose short intense pulse with £2 » £2[, ~ F, 
which is used for population transfer. 
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FIG. 3: (Color online) Maximum excitation probability / > ^ Iax as a 
function of the mean photon number N with the initial coherent state 
Gaussian pulse for different bandwidths. C1' = 2AT is the optimized 
bandwidth for a Gaussian pulse (c.f. Fig[2]l. 



B. Pulse shape effects 



In general, the excitation probability depends on the 
specific shape of the input pulse. Here we studied the 
following six pulse shapes, see Table U 



TABLE I: Definition of pulse shapes 



Type of pulse 



Wave function for pulse bandwidth 



Gaussian pulse 
Hyperbolic secant pulse 

Rectangular pulse 

Symmetric exponential pulse 



f(f) = (g) exp(-^ 2 ) 
JfsechCOf) 



forO < t < 



Decaying exponential pulse £(f) : 



Rising exponential pulse 



I 0, else 

f(f) = Vfiexp(-n|f|) 

Vnexp(-fr), for/>0 
0, forf<0 

Vfiexp(ff), for/<0 
0, forf>0 



For single photon Fock state, the excitation probability 
has a peak value of about 0.8 with optimum bandwidth for 
the first four pulse shapes, shown in Fig|4](a)-(d), indicating 
that the photon absorption is less sensitive to pulse shape 
effects such as discontinuities. For the decaying exponential 
pulse, the maximum excitation probability is only 0.54, see 
FigH](e). A particularly interesting case may be that of the 
rising exponential single photon Fock state pulse, shown in 
FigS (f), for which the corresponding maximal excitation 
probability is 0.995 with a optimal bandwidth of Qo = IF. 
This agrees well with the prediction that for the aim of unit 
absorption probability, the incident photon must possess the 
time reversed properties of the spontaneously emitted photon. 
Since the spontaneous decay is exponential, the temporal 
envelope of the pulse has to be rising exponential lfl5l[l8ll . 

On the other hand, for an initial single photon coherent 
state pulse with optimum bandwidth, the maximum excitation 
probability is much lower, around 0.48 for the first four 
pulse shapes and 0.4 and 0.56 for the decaying and rising 
exponential pulse, respectively. Apparently the excitation is 
more efficient if exactly one photon is present instead of a 
distribution with mean one. This emphasis the importance of 
generating single photon source rather than using attenuated 
laser pulse in applications where a high absorption is desired. 

For the explicit values of optimum bandwidth needed to 
achieve maximum excitation probability, see Table. |ll] 
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FIG. 4: (Color online) Excitation probability P e (t) as a function of 
time for the coupling strength g{t) given in Eq.J22t with r p = T (A = 
87r/3). The single photon Fock state pulse with optimal bandwidth is 
shown in grey; the corresponding excitation probability is given by 
the solid black line. The dashed blue line represents the excitation 
probability for a single photon coherent state pulse of a similar 
shape but different (optimized) bandwidth, (a) Gaussian pulse, 
(b) Hyperbolic secant pulse, (c) Rectangular pulse, (d) Symmetric 
exponential pulse, (e) Decaying exponential pulse and (f) Rising 
exponential pulse. 



TABLE II: Optimum bandwidth and maximum excitation 
probability, the results (*) were also obtained in Ref. lT^I with 
a different method. 



Tape of pulse 


State 


Optimum Q./T 


Maximum P e (t) 


Gaussian pulse 


Iff) 


2.4 


0.48 


U> 


1.5 


0.80 * 


Hyperbolic secant pulse 


Iff) 


2.0 


0.48 


ID 


1.3 


0.80 


Rectangular pulse 


Iff) 


1.3 


0.48 


U> 


0.8 


0.81 


Symmetric exponential pulse 


Iff) 


1.4 


0.48 


ID 


0.9 


0.79 


Decaying exponential pulse 


Iff) 


1.4 


0.37 


ID 


1.0 


0.54 * 


Rising exponential pulse 


Iff) 


1.9 


0.56 


ID 


1.0 


0.995 * 



C. Damped Rabi oscillation 

In Fig|5] the probability of exciting the atom for an initial 
coherent state Gaussian pulse is evaluated for various mean 



photon numbers N = (1,10,50). For large mean photon 
number, damped Rabi oscillations are observed. In the limit 
of very large mean photon number, one would recover the 
textbook predictions forclassical light pulses" |p23l p. 151]. 
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FIG. 5: (Color online) Excitation probability P e (t) as a function of 
time for initial coherent state Gaussian pulse with optimal bandwidth 
Q.' = 2.4T for different mean photon numbers N. 



V. DISCUSSION OF REALISTIC FOCUSING 

Finally we present a brief review of ongoing experiments 
in order to consider the excitation probability in realistic tight 
focusing configurations. 

In the case of a parabolic mirror with a half opening 
angle of 134° as it is used in the experiment described in 
Refs. lfl3l[l4ll . the corresponding weighted solid angle reaches 
A = 0.94 x 87r/3, and thus one may achieve a maximal 
excitation probability of 0.94 with rising exponential shape 
for a single photon Fock state pulse, 0.54 for a single photon 
coherent state pulse, 0.75 for a Gaussian single photon Fock 
state pulse and 0.46 for a single photon coherent state pulse. 

In Ref.ll, Gil, a high aperture lens with NA = 0.55 and 
/ = 4.5 mm is used to focus down a Gaussian beam. The 
weighted solid angle depends on the focusing strength u := 
wi/f, where is the beam waist. A maximum overlap 
of A = 0.364 x 8^/3 is expected at focusing strength u = 
2.239. With a rising exponential shape, we predict a maximal 
excitation probability of 0.36 for a single photon Fock state 
pulse and 0.27 for a single photon coherent state pulse. For a 
Gaussian shape, we predict a maximal excitation probability 
of 0.29 for a single photon Fock state pulse and 0.23 for a 
single photon coherent state pulse. 



VI. CONCLUSION 

In conclusion, with the help of time dependent 
Heisenberg-Langevin equations, we studied the interaction 
between a single two-level atom and a propagating pulse at 
the quantum level. We have presented a general approach 
and a scalar model to treat the atom-pulse interaction. For 
strong focusing configurations we account for the overlap of 
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the incoming pulse mode with the respective atomic dipole 
pattern. 

The effect of temporal-spectral features of the single photon 
Fock state pulse and coherent state pulse on excitation 
probability of the atom has been investigated. With Gaussian, 
hyperbolic secant, rectangular and symmetric exponential 
shape pulses, the achievable maximum excitation probability 
is around 0.8 for single photon Fock state pulse and 0.48 for 
single photon coherent state pulse. More importantly, with a 
rising exponential shape, the maximum value is nearly 1 for a 
single photon Fock state pulse and 0.55 for a single photon 
coherent state pulse, which is in agreement with the time 
reversal argument. As an example, the effect of bandwidths 
and mean photon numbers of a Gaussian pulse is analyzed. 
We also survey some current arrangements designed to couple 



photons and atoms in order to assess their potential for high 
excitation probabilities. The presented model provides a 
suitable foundation to further study the pulse shape effects 
of quantized continuous light fields in both scalar and full 
three-dimensional treatment. 
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